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Within the Kubo-Greenwood formalism we use the fully relativistic, spin-polarized, screened 
Korringa-Kohn-Rostoker method together with the coherent-potential approximation for layered 
systems to calculate the resistivity for the permalloy series NicFei_c. We are able to reproduce 
the variation of the resistivity across the entire series; notably the discontinuous behavior in the 
vicinity of the structural phase transition from bcc to fee. The absolute values for the resistivity 
are within a factor of two of the experimental data. Also the giant magnetoresistance of a series of 
permalloy-based spin- valve structures is estimated; we are able to reproduce the trends and values 
observed on prototypical spin- valve structures. 

PACS numbers; 75.30.Gw, 75.70.Ak, 75.70.Cn 
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I. INTRODUCTION 

Permalloy is perhaps the most commonly used mag- 
netic material in a variety of devices. This arises primar- 
ily from its magnetic properties which combine a rea- 
sonably high magnetic moment with a low coercivity. In 
many of its applications one passes a current through the 
material, so that its electrical transport properties are 
also important. There have been several calculations of 
the transport properties of permalloy, but none has been 
completely successful. Therefore, one of the outstand- 
ing challenges has been an ah initio calculation of the 
electrical resistivity and anisotropic magnetoresistance 
across the entire NicFei_c series. Here we present the 
first part which addresses the resistivity of bulk permal- 
loy and also the giant magnetoresistance of spin-valve 
structures containing layers of permalloy. As the lattice 
structure changes from bcc to fee at c ~ 35% we present 
results for both lattice structures. We are able to re- 
produce the overall variation of the resistivity across the 
series; in particular we find the signature changes in re- 
sistivity about c ~ 35% as the lattice structure changes 
from bcc to fee. Remarkably we find the magnitude of 
the resistivity is within a factor of two of the experimen- 
tal values; noteworthy our calculated values are larger, 
so that when one relaxes the approximations we made 
one reduces the resistivities and approaches the measured 
ones. 

Here we outline our ah initio calculation and present 
our results for the resistivity across the entire NicFei_c 



series. Then we conclude with results on prototyp- 
ical spin-valve structures containing permalloy layers, 
e.g., NicFei_c(100A) / Co(6A) / Cu(9A) / Co(6A) / 
NicFei_c(5A). Here we find current in the plane of the 
layers (CIP) magnetoresistance ratios that compare fa- 
vorably with those observed in these structures. 



II. METHOD OF CALCULATION 

In terms of the Kubo-Greenwood approach for disor- 
dered layered systems with growth direction along the 
sucface normal (2;-axis) the in-plane conductivity is given 

n 

cr^^(n;c;C) = ^ a'^^{c;C), fie {x,y}, 

(1) 



ij _ h 



Tr(j^ImG+(eF)J;^ImG+(eF)) 



where n denotes the number of layers considered, c = 
{ci, C2, . . . , c„} is a set containing the layer- wise compo- 
sitions, e.g., the layer- wise concentrations in an inhomo- 
geneously disordered binary alloy system Ac-Bi-c^ (in a 
homogeneously disordered alloy Ci = c, V«; whereas in a 
spin- valve system A^ / B„ / A^ the Ci = 1, Vi G A^, and 
Ci = 0,Vz G B„). In principle the magnetization in each 
layer of a multilayered structure can be different, and 
C ~ {Ml, M2, . . . , M„} denotes a particular magnetic 
configuration, where the M; refer to the orientations of 
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the magnetizations in the individual layers, fi^/^ is the 
conductivity that describes the current in layer i caused 
by an electric field in layer j, A'o the number of atoms per 
plane of atoms, fiat is the atomic volume, (• • •) denotes an 
average over configurations, is the fi-th component of 
the current operator referenced to the i-th plane, and 
is the electron propagator (one-particle Green's function) 
from plane i to j at the Fermi energy ep. The resistiv- 
ity corresponding to the conductivity given by Eq. (|l|) is 
then defined by 

Pf,f,{n;c;C) = l/(Tf,f,{n;c;C), e {x,y}. (2) 

By assuming two different magnetic configurations, C and 
C, the relative change in the resistivities between them 
is given by the ratio 

p^^(n;c;C') 

For example in spin-valve structures C and C refer to 
"antiparallel" and "parallel" alignments of the magneti- 
zations of the magnetic slabs, and this ratio is referred to 
as the giant magnetoresistance (GMR) ratio. The ratio 
in Eq. ^ also serves as a definition for the anisotropic 
magnetoresistance of homogeneously disordered bulk al- 
loys, Ci = c, Vi, with uniform magnetization, M.; = M, Vi, 
with C and C referring to configurations with the magne- 
tization pointing uniformally parallel and perpendicular 
to the current axis (however, in this case the ratio is usu- 
ally taken with respect to an averaged resistivity) . 

In the present paper all calculations are based on 
self-consistent effective potentials and effective exchange 
fields as obtained previouslynl by using the fully relativis- 
tic, spin-polarized, screened Korringa-Kohn-Rostoker 
method together with the coherent-potential approxima- 
tion for layered systems. For permalloy a fully rela- 
tivistic spin-polarized calculation of the resistivity is es- 
sential in order to properly account for the spin-orbit 
coupling of the electron states as well as for the scat- 
tering|-.bjy the disordered constituents of these binary 
alloys.0cl The surface-Brillouin-zone integrals needed in 
the evaluation of the electrical conductivity within the 
Kubo-Greenwood approacnJ were obtained by consider- 
ing 1830 k|| -points in the irreducible wedge of the surface 
Brillouin zone. All scattering channels up to and includ- 
ing £niax = 2 were taken into account. Here we have cho- 
sen the magnetization to be uniformly perpendicular to 
the planes of atoms, i.e., the magnetization points along 
the 2;-axis. With this choice the xx- and yy-components 
in the above equations are identical, hence we have taken 
the direction of the current (electric field) to be along the 
cc-axis so as to simulate the CIP transport geometry (with 
the magnetization pointing uniformally perpendicular to 
the current axis). 

Consider the following systems: 

(a) Ni,Fei_e(100) / (Ni,Fei_e)„ / Ni,Fei_,(100) 

(b) Ni,Fei_c(100) / (Ni,Fei_c)„ / Vacuum (4) 

(c) Vacuum / (NicFei_c)„ / Vacuum . 



All three cases refer to the in-plane conductivity for a 
film with n monolayers of permalloy; what distinguishes 
them arc the boundary conditions.El In case (a) we have 
outgoing boundary conditions on both sides as there is 
nothing that separates the film from the substrate it is 
deposited on nor from the overlayer that caps it, i.e., 
some electrons will leak out into the semi-infinite leads 
and their momentum information will be lost hence pro- 
ducing an additional contribution to the resistivity. In 
(b) we have a free surface on one side which makes it 
refiecting (the workfunction is finite so that the potential 
barrier is not infinitely high and there will be a slight 
leaking out of electrons into the vacuum that has to be 
taken into account). In the third case (c) we have a 
perfectly flat free standing film with reflecting boundary 
conditions on both sides. Only in (c) the calculated in- 
plane conductivities/resistivities for a perfectly flat film 
are principally independent of the film thickness n, at 
least for n > 12; roughness would act like an outgoing 
boundary condition. However, for all three cases one 
obtains the so-called bulk residual resistivity p°{c;C) of 
substitutionally disordered binary alloys with respect to a 
particular magnetic configuration by taking the infinite- 
thickness limitQ 

p"(c;C)= lim /9^^(n;c;C). (5) 

n — >oo 

This infinite-thickness limit (n —>■ oo) can be obtained 
easily since np^^(n; c; C) becomes linear in n, if n is large 
enough; its slope simply refers to the extrapolated resis- 
tivity (see also the discussion of Fig. 2 in the next Sec- 
tion). 

In addition, for computational purposes a finite imag- 
inary part 5 to the Fermi energy has to be used in- the 
calculation of the conductivity (J^^{n\ c;C;eF + i6)B To 
designate resistivities calculated with finite S we adopt 
the notation 

Pfj,fj,{n;c;C;S) = l/CT^^(n; c;C; ef + «<5), p. e {x,y}- 

(6) 

The actual resisitivity is arrived at by numerically taking 
the limit as the imaginary part goes to zero 

Pf,f,{n;c;C) = lim p^^(7i; c; C; (5). (7) 

If we are interested in the bulk residual resistivity the 
{S — * 0) extrapolation is simplified by the fact that the 
slope of npfj_fj_{n]c;C;S) (determined at large enough n) 
is itself linear in 6. 

As a check on the above extrapolation procedure we 
used to arrive at bulk values of the resistivity from those 
on films of finite thickness (Eq. (||)) we have calculated in- 
plane resistivities (with the magnetization pointing uni- 
formly along the z-axis) for 30 monolayers of permalloy 
with a homogeneous Ni concentration of c = 0.85 for a 
fixed 6 = 2 mRyd with different boundary conditions: 
For the reflecting boundary conditions (case (c)) we find 
a resistivity of p'^^{30; 0.85; z; 2) = 20.6 pQcm. A similar 
calculation with outgoing boundary conditions (case (a)) 
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yields p°;i*(30; 0.85; z; 2) = 31.5 nflcm which by taking 
n ^ oo yields p°"'(cxd; 0.85; z; 2) = 20.6 /xficm; precisely 
the same value as found for the reflecting boundary con- 
ditions (case (c)). This unequivocally demonstrates that 
the procedure outlined above is able to remove the resis- 
tivity that_arises from the boundary conditions on finite 
structures .0 In Fig. 1 we show as an illustration the vari- 
ation of the layer-wise conductivities 

n 

i=i 

across the film with reflecting and outgoing boundary 
conditions. For the outgoing boundary conditions where 
some electrons leak into the leads all layer-wise conduc- 
tivity contributions are substantially smaller than for 
the reflecting boundary conditions, the difference becom- 
ing progressively larger when we get closer to the outer 
edges of the film. When we apply reflecting bound- 
ary conditions the layer-wise conductivities are practi- 
cally constant with small oscillations from the bound- 
ary towards the center of the film resembling the well- 
known Friedel oscillations. When we take the result 
for (5 = 2 mRyd corresponding to outgoing boundary 
conditions to — > oo (Eq. (p|)), and then numerically 
take the limit as (5 (Eq. (0)), we find p°(0.85;z) = 
/9™*(oo; 0.85; z; 0) = 7.1 fiilcm which is the calculated 
resistivity for Nio.85Feo.15 to be compared with experi- 
ment. It should be noted, however, that the extrapo- 
lation procedures give an upper limit for the resistivity 
since due to computational constraints our calculations 
can be done only up to n ~ 45 with S > 2 mRyd (i.e., 
the slopes of Tip^^{n;c;C;S) might be slightly too high, 
and close to (5 = there might be small deviations from 
linearity). Parenthetically, a similar extrapolation pro- 
cedure on a film with no disorder, i.e., c — 1.0, yields 
finite Pxx{n] 1.0; C; 5) in both cases (a) and (c) for finite 
5; however, the (n ^ 00, (5 ^ 0) extrapolated value 
p°"*(cx); 1.0; C; 0) for outgoing boundary conditions, as 
well as the [5 0) Hmit p^'?^(n; 1.0; C; 0) for reflecting 
boundary conditions, are both zero. 

III. RESULTS 

A. Bulk alloys 

For case (b) of Eq. (^, Eq. (^ is precisely the pro- 
cedure pursued in experimental studies; see Fig. 6 of 
McGuire and PotterQ where the thickness dependence of 
the resistivity of Nio.82Feo.i8 fllms is shown. It is well 
known that the film resistivity increases as films become 
thinner; the resistivity noticeable increases when the film 
thickness becomes much smaller than the mean free path 
of the conduction electrons. It is quite interesting to re- 
alize that the measured resistivity reaches its asymptotic 
bulk value only for film thicknesses well above 500A. In 
Fig. 2 we display the thickness dependence of our calcu- 
lated p°"*(7i; c; z; 0) for case (a) of Eq. (^) as a function 



of the number of layers n considered in the summation 
in Eq. (|l|). As can be seen, for small values of n the cor- 
responding changes in /9™*(n; c; z; 0) are large, whereas 
— ^just as for the experimental data — for large enough n 
/9°"*(n; c; z; 0) approaches towards its asymptotic value, 
i.e., the functional behavior of n/9°"*(n; c; z; 0) becomes 
linear in n. Inspecting our (n —>■ 00) extrapolation we 
find (again in accordance with experiment) that above 
585A (= 330 monolayers) the deviations from the bulk 
residual resistivity are smaller than 1 fiflcm, i.e., in or- 
der to talk about "bulk" in terms of resistivities one has 
to consider a tremendously thick film of several thou- 
sands A. It should be noted that in order to arrive at 
bulk resistivities the same kind of linear extrapolation 
procedure to very large n (Eq. (|^)) is done both with 
experimental data and with our calculated resistivities 
(Eqs. (0) and (|)). 

In Fig. 3 we show our theoretical asymptotic bulk resis- 
tivities p°(c; z) = /o°"*(oo; c; z; 0) obtained by performing 
both limits (n —^ 00 and (5 — > 0) for the entire concentra- 
tion range tog^her with low temperature (4.2 K) exper- 
imental dataElEl for Ni concentrations above 70%. The 
overall trends of measured data can be obtained from an 
inspection of Fig. 1 of the section on NicFei_c in Ref. |l^ 
which shows the temperature dependence of the resis- 
tivity for the whole range of concentrations. From this 
figure one can see that by lowering the temperature the 
electrical resistivity becomes increasingly sensitive to the 
structural phase transition from bcc to fee; at the temper- 
ature of liquid nitrogen (—195 C, the lowest isothermal 
curve displayed) in the Fe-rich bcc a-phase the resistivity 
reaches a maximum at about 15% of Ni and then slowly 
decreases up to 30% of Ni, whereas in the Ni-rich fee 7- 
phase the resistivity starts to grow below 50% of Ni and 
seems to diverge near the critical concentration for the 
structural phase transition. In the —195 C curve we also 
discern a weak shoulder at about 85% Ni; this feature is 
well regiolved in the low temperature (4.2 K) experimental 
data.Lru As can be seen from Fig. 3 our theoretical values 
reproduce very well the vicinity of the structural phase 
transition from bcc to fee at about 35% of Ni. There is 
indeed a maximum in the bcc phase at about 15% of Ni; 
the onset of the phase transition is clearly visible. Also 
— ^just as the experimental data — in the fee phase the 
theoretical values show a kind of shoulder at about 85% 
of Ni. 

For this concentration (c = 0.85) we performed an ad- 
ditional calculation using a maximum angular momen- 
tum quantum number £niax — 3 (see also Ref. |l^) . As can 
be seen from Fig. 3, the difference between the imax = 3 
and the iniax — 2 (that is used in all other calculations) 
results is rather small; however the trend for the inclu- 
sion of higher scattering channels is to lower (slightly) 
the resistivity. 

Although the theoretical resistivities agree rather well 
with the experimental ones, the difference at a particu- 
lar concentration being of the order of 3 /if2cm, it seems 
appropriate to comment on the putative sources of the 
difference. First of all, in the present use of the Kubo- 



Greenwood equation for layered systems no vertex cor- 
rections arising from the configurational average of the 
product of two Green's functions arc included isee also 
the discussion in Ref . |l|) . In another calculationO'cl of the 
resistivity of bulk NicFei_c the vertex corrections were 
calculated and found to—jpe quite small over the entire 
range of concentrationSjEj i.e., no larger than the differ- 
ences between the ^max = 2 and fmax = 3 results at 
c = 0.85 in Fig. 3. Second, the single site approxima- 
tion to the coherent-potential approximation is used to 
describe the electronic structure of substitutional alloys. 
This in turn implies that short range order or concen- 
tration fluctuations are excluded. Both the inclusion of 
vertex corrections and the inclusion of short range order 
will reduce the present theoretical resistivity values. Fur- 
thermore, in magnetic alloys magnetic ordering or clus- 
tering can occur; as it implies some ordering, which we 
have not taken into account, it will again (slightly) lower 
the resistivity. |-||-| 

In two previous theoretical papergaa the validity of 
the so-called two-current model for transport in fer- 
romagnetic systems with strong spin-dependent disor- 
der was tested; the resistivity and the anisotropic mag- 
netoresistance for NicFei_c bulk alloys was calculated 
by solving the Kubo-Greenwood equation using a spin- 
polarized, relativistic version of the Korringa-Kohn- 
Rostoker method together with the coherent-potential 
approximation. As these calculations were done on a 
fee lattice for the whole concentration range their results 
are not applicable to the Fe-rich bcc alloys. Nonetheless 
a very weak shouJder at about 80% Ni can be discerned 
from their Fig. 413 — in agreement with the experimental 
data and our present calculation, see Fig. 3. However, 
the resistivities they find for 80% Ni, p ~ 1.8 fiftcm, as 
well as the maximum at about 40% Ni, p ~ 3.2 /ii7cm, are 
too low by about a factor of two as compared to available 
experimental data. 

B. Spin- valve structures 

The main purpose of our study has been to show how 
well resistivities for disordered magnetic systems with 
only two-dimensional translational symmetry (films) can 
be described, and to confirm that by taking the appro- 
priate limit we derive the bulk resistivities from those 
calculated for films. Since a large number of spin- valve 
structures that are used in applications, e.g., the read 
heads of hard disk drives, contain permalloy layers, we 
conclude with some results on their transport proper- 
ties. Including permalloy layers into a multilayer sys- 
tem the electronic structure and the magnetic proper- 
ties are considerably more complicated than for bulk 
permalloy; specifically they are sensitive to the thick- 
ness and Ni concentration c of the permalloy layers, and 
to the constituency of the neighboring layers. There- 
fore for concreteness we focus on a permalloy-based 
system that has been studied, e.g.. Si-substrate / Cu- 
contacts(1500A) / NicFei_c(100A) / Co(6A) / Cu(9A) / 



Co(6A) / NicFei_c(5A) / NiO(250A),El and we replace it 
with the reference system NicFei_c(100) / (NicFei_c)i2 
C04 Cu„ C04 (NicFei_c)3 / Vacuum. In all cases shown 
a fcc(lOO) stacking was chosen, the lattice spacing refers 
to that of the corresponding fee bulk NicFei_c system. 
In this transcription the rather thick permalloy layer of 
about lOOA (= 56 monolayers) is considered as a sub- 
strate, the insulating NiO layer is replaced by a free sur- 
face (vacuum), and all other thicknesses are expressed (in 
monolayers) as close as possible to the measured struc- 
ture. The 12 layers of NicFei_c are deliberately chosen 
so as to guarantee a smooth matching to the substrate 
via a selfconsistent calculation; this in turn would also 
allow us to change the left boundary condition such that 
in principle a system of the type Vacuum / (NicFei_c)5g 
C04 Cu„ C04 (NicFei_c)3 / Vacuum could be described 
by assuming that the layers in the very interior of the 
thick NicFei_c layer exhibit bulk-like properties. 

We have considered two series: in the first the Cu- 
spacer thickness is confined to n = 5 and we vary the Ni 
concentration c while maintaining an fee lattice through- 
out; in the second we set c = 0.85 and vary the number of 
Cu-spacer layers n. While the structure of NicFei_c bulk 
alloys is bcc for c < 35% when there are relatively thin 
layers in a fee multilayered structure it is appropriate to 
maintain the fee lattice structure even for the Fe-rich al- 
loys. In the parallel configuration of the spin valve the 
magnetizations of all permalloy and Co layers are aligned 
uniformly perpendicular to the plane of the layers; in 
the antiparallel configuration they are still normal to the 
layers but the magnetizations of the two magnetic slabs 
(separated by the Cu-spacer) point in opposite directions 
to one another. The conductivity is evaluated at a finite 
imaginary part to the Fermi energy, S — 2 mRyd; this 
is necessary for otherwise the calculations would take an 
inordinate time. In principle, in order to obtain the re- 
sistivity we have to take the limit as (5 — > (see Eq. (|^)). 
However, it is virtually impossible for technical reasons 
to do this up till now for a whole series of spin-valve 
structures; therefore we do not show the resistivities we 
calculated for finite i5. But with the assumption that 
both resistivities for the parallel and the antiparallel con- 
figuration of the spin valve scale uniformly with S, and 
therefore cancel out in the ratio R (Eq. (||)), we have 
used them to estimate the CIP-MR ratios that one can 
expect from these spin- valve structures. 

In Fig. 4a we show the variation of the calculated CIP- 
MR ratio of the reference system with the Ni concentra- 
tion c for the fixed Cu-spacer thickness n — 5. As can be 
seen the CIP-MR varies only very little (between about 
17.5-21.5%) over the whole range of concentrations con- 
sidered. In Fig. 4b we show the variation of the CIP-MR 
ratio for the reference system corresponding to c = 0.85 
with respect to the number of Cu-spacer layers n. By 
doubling the number of Cu-spacer layers, i.e., by going 
from n = 4 to n = 8, the CIP-MR is reduced by almost 
a factor of two. This dramatic decrease can be under- 
stood in terms of a simple model: since on increasing the 
thickness of the nonmagnetic spacer layer tnm the CIP- 
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MR decreases as (1/inm) exp(— inm/Anm),t£l where Anm is 
the mean free path of the conduction electrons in the non- 
magnetic spacer layer, the ratio Rxxit^m) / Rxx{'2tnm) ~ 2 
for tnm ^ Anm (which is a necessary condition to find 
at all a CIP-MR and certainly valid for our tnm vary- 
ing between 7.1 A and 14.2 A). Our findings seem tq 
be confirmed by experiments made on similar systemsEj 
albeit at room temperature, while our calculated values 
are for T = K. Depending on the actual thicknesses 
of: the bottom and top NicFei_c layers, of the Co slabs 
and of the Cu spacer, and the pinning and capping lay- 
ers used, the experimental CIP-MR ratios in permalloy- 
related spin-valve systems have been found to be in the 
range between 10 and 20%. Fig. 4 nicely shows that this 
is about the range that our theoretical calculations pre- 
dict when we use the resistivities calculated with finite 6, 
and when we make the assumption that for both the par- 
allel and the antiparallel configurations of the spin valve 
the resistivity scales uniformly with S and therefore can- 
cels out in the CIP-MR ratio R. 



IV. SUMMARY 

By using our fully relativistic, spin-polarized, screened 
Korringa-Kohn-Rostoker method together with the 
coherent-potential approximation for layered systems we 
have calculated the resistivity across the entire NicFei_c 
(permalloy) series of bulk alloys within the Kubo- 
Greenwood approach. We are able to reproduce the over- 
all variation of the resistivity as a function of the Ni con- 
centration; in particular we reproduce the discontinuity 
in the vicinity of the structural change from bcc to fee. 
Also we obtain very reasonable absolute values for the 
resistivities. We have also estimated the CIP-MR ratios 
for permalloy-based spin- valve structures and within the 
limits of our current calculation find ratios in reason- 
ably good agreement with those measured on prototypi- 
cal spin valves; there is of course the caveat that most of 
the data on spin valves is quoted at room temperature 
while our results are for T = K. It should be stressed 
that our results are based on ab initio calculations with 
no adjustable parameters. Viewed in this light the cal- 
culated resistivities and GMR ratios are noteworthy. 
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FIG. 1. Variation of the layer-wise conductivities 
o"a;a;(30; 0.85; z; tF + *2) (Eq. (^) across the film for 30 mono- 
layers of permalloy with a homogeneous Ni concentration of 
c = 0.85 for a fixed 5 = 2 mRyd with different boundary 
conditions: The values for the reflecting boundary conditions 
are show as empty squares, those for the outgoing boundary 
conditions as full circles. 



FIG. 3. Concentration dependence of the current-in-plane 
resistivity of "bulk" NicFei_c alloys calculated as the in- 
finite-thickness limit of NicFei-c(lOO) / (NicFei_c)„^^ / 
NicFei_c(100) (p°(c; z) = p°t(.oo; c; z; 0)). The results for the 
bcc Q-phase are shown as empty squares, for the fee 7-phase 
as full circles, the open circle at c = 85% refers to a calcula- 
tion with linsx = 3. Low temperature (4.2 K) experimental 
values displayed by full diamonds are taken from Refs. 
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FIG. 2. Thickness dependence of the current-in-plane re- 
sistivity in bcc (a) and fee (b) NicFei_c(100) / (NicFei_c)„ 
/ NicFei_c(100) alloys (p°"*(n; c; z; 0)). The Ni concentration 
in % is indicated explicitly, n denotes the number of layers 
considered in the summation in Eq. mh. 
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FIG. 4. Current-in-plane magnetoresistance for permal- 
loy-based spin valves: (a) Concentration dependence for 
fee-based NicFei_c(100) / (NicFei-c)i2 C04 Cus C04 
(NicFei_c)3 / Vacuum systems, (b) Dependence on the num- 
ber of Cu-spacer layers n in fcc-based Nio.g5Feo.i5(100) / 
(Nio.ssFeo. 15)12 Cu„ C04 (Nio.85Feo.i5)3 / VacuuirkT-|E«- 



perimental values are found to be in the range 10-20%.ll3i 



